Critical behavior of two freely evolving granular gases separated 

by an adiabatic piston 

J. Javier Brey and Nagi Khalil 
Fisica Teorica, Universidad de Sevilla, 
Apartado de Correos 1065, E-^1080, Sevilla, Spain 
(Dated: September 24, 2010) 

Abstract 

Two granular gases separated by an adiabatic piston and initially in the same macroscopic state 
are considered. It is found that a phase transition with an spontaneous symmetry breaking occurs. 
When the mass of the piston is increased beyond a critical value, the piston moves to a stationary 
position different from the middle of the system. The transition is accurately described by a simple 
kinetic model that takes into account the velocity fluctuations of the piston. Interestingly, the final 
state is not characterized by the equality of the temperatures of the subsystems but by the cooling 
rates being the same. Some relevant consequences of this feature are discussed. 

PACS numbers: 45.70.-n,05.20.Dd,5.40.-a 
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I. INTRODUCTION 

The so-called adiabatic piston []] is a typical example showing the relevance of fluctuations 
to explain the macroscopic behavior of a system. It consists of an isolated cylinder containing 
two gases separated by an adiabatic piston (no internal degrees of freedom). Usually, the 
system is initially prepared with the gases in both compartments at independent equilibrium 
states and the piston fixed by a clamp at a given position. Then the clamp is removed and 
the piston is free to move without friction with the container. One of the interesting features 
of this system is that equilibrium thermodynamics cannot predict the final position of the 
piston and the states of the gases to both sides of the piston. Nevertheless, when fluctuations 
are taken into account, it follows that the piston moves until the system relaxes to mechanical 
and thermal equilibrium with equal pressures and temperatures in both compartments. In 
the last years, the problem has attracted a lot of attention 2j-|4j, mainly stimulated by the 
seminal paper by Lieb j^j and the suitability of the model to investigate fundamental issues 
in mesoscopic systems. An illuminating review of the adiabatic piston is given in ref. 6]. 

Granular gases, modeled as ensembles of particles colliding inelastically, exhibit many 
similarities with molecular gases, but also a rather large number of peculiar behaviors as a 
consequence of energy dissipation. In particular, they do not have equilibrium states and 
are inherently non-equilibrium systems. For this reason, they have been considered as a 
proving ground for kinetic theory and non-equilibrium statistical mechanics [?], Q • In this 
context, an adiabatic piston separating two granular gases appears as a natural system to 
investigate the effects of the irreversible dynamics on the interaction between macroscopic 
systems. Brito et al. j^] considered the case in which the two inelastic gases are initially 
prepared in the same macroscopic state, the homogeneous cooling state, and the particles 
collide elastically with the piston. They found that the piston eventually collapsed to one 
of the sides of the container. This behavior was explained in terms of the instability of the 
initial state, leading to an spontaneous left-right symmetry breaking. 

Asymmetric inelastic pistons, with both sides made of different materials with different 
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inelasticities, have also been considered [10|, LL3J. In these studies, it has been assumed 
that the piston moves in an infinite bath of elastic particles either with a Gaussian velocity 
distribution 1 10] or with a velocity distribution in which only two discrete opposite velocities 
are possible Of course, the two kind of systems, elastic piston and inelastic gases versus 
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infinite baths of elastic gases colliding inelastically with the piston, are rather different. In 
the former case, the system is continuously cooling while in the latter a steady state is 
reached after some transient. 

Here the model of the symmetric elastic piston separating two inelastic gases is revisited. 
The aim is to incorporate in the description given in |9[ the effect of the velocity fluctuations 
of the piston. This is the mechanism for which energy is transferred between the left and right 
hand side gases, leading in the case of molecular gases to an equilibration of the temperature. 
It will be shown that some kind of equilibration also happens in the case of granular gases, 
although driven by a parameter different from the temperature. The quantity that becomes 
the same is the cooling rate of the two gases and also of the piston. Quite interestingly, these 
conditions lead to the existence of a non-equilibrium phase transition with an spontaneous 
symmetry breaking, when the mass of the piston exceeds some critical value, which depends 
on the number of particles and on the inelasticity of the collisions between them. The phase 
transition is accurately described by a simple kinetic theory that incorporates the effect of 
the velocity fluctuations of the piston. The results fit into an scenario in which the whole 
system composed by the gases in both compartments and the piston reaches a cooling state 
characterized by a unique temperature parameter. The results of ref. (9] are recovered in 
the limit of an infinitely massive piston. 

The remaining of the paper is organized as follows. In the next section the system under 
consideration is specified and the simple theoretical model is formulated. Also the states 
to be addressed are specified. Evolution equations for the temperature parameters of the 
granular gases and the piston in those states are written down. Two different kind of states 
are possible. In one group, the system is symmetric to both sides of the piston, while in 
the other group the initial left-right symmetry is spontaneously broken. Both states are 
discussed in Sees. II I II and IIV} respectively. A control parameter characterizing whether the 
state is symmetric or asymmetric is introduced. In Sec. |V]the theoretical predictions are 
compared with molecular dynamics simulations and a good agreement is observed in all the 
accesible parameter region. The final section contains some concluding remarks and also a 
discussion of the possible implications of the reported results for the understanding of the 
interactions between macroscopic granular systems. 
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FIG. 1: Sketch of the adiabatic piston considered here. Each of the compartments contains iV 
inelastic hard spheres or disks of mass m and diameter a. The piston can move without friction 
along the X clXIS, clS indicated in the figure. 

II. THE MODEL 

The system considered consists of a cylinder of length L x filled with 2N inelastic hard 
spheres (d = 3) or disks (d = 2) of mass m and diameter a. The cylinder is divided into two 
compartments by a piston of mass M perpendicular to the symmetry axis of the container, 
taken as the x-axis (Fig. [T|). The piston can move without friction along the x direction, its 
motion being entirely induced by the collisions with the gas particles. There are the same 
number of particles, N, at each side of the piston. 

The system is isolated and evolves freely in time. Collisions between particles are char- 
acterized by a constant (velocity-independent) coefficient of normal restitution a, defined in 
the interval < a < 1. On the other hand, collisions of the particles with the piston and 
with the walls of the container are elastic. Attention will focus here on states in which the 
position of the piston has a stationary value xp. In these states, mechanical equilibrium 
requires that the hydrodynamic pressure p be the same at both sides of the piston. The vari- 
ables corresponding to each of the compartments will be identified by a subscript % = 1,2, 
respectively. It will be assumed that the gas remains always very dilute in both compart- 
ments and that, as a lowest order approximation, the granular gas is homogeneous inside 
each compartment and it can be described by the equations for the homogeneous cooling 
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state (HCS) Q, [ijj]. This implies that the relative flux of energy through the piston is very 
low as compared with the characteristic relaxation time of the HCS. Also, it requires that 
the system is not too large so that the HCS is stable 13j and the shearing and clustering 
instabilities are avoided. Then, for the states under consideration, it is Pi(t) = P2{t) or, 
equivalently, 

n l T 1 (t)=n 2 T 2 (t), (1) 

where n and T denote the number density and the granular temperature respectively. Note 
that while the temperature and the pressure will be time dependent, the stationarity of xp 
implies the same property for the density, since 

N N 



n 2 



(2) 



Sx P S(L x -x P ) 
Here S is the section (area for d = 3 and length for d — 2) of the piston and the container. 



Energy balance equations for each of the compartments when the piston is at rest are given 
by 

J rIT- J 

(3) 



< L N d Il = - d L N C !l T i + Q l S 1 
2 dt 2 s ^ 



i = 1,2. The first term on the right hand side is the Haff 's law describing the cooling of the 
iCS due to the inelasticity of collisions The explicit expression of the cooling rate is 
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with rjo being the elastic (a = 1) shear viscosity, 
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and C*(°0 a dimensionless function of the restitution coefficient, 
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The above expressions 
Sonine approximation 



9 + 2Ad + (8d - 41)a + 30a 2 (l - a) ' 
rave been derived from the Boltzmann equation in the so-called first 
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The second term on the right hand side of Eq. ([3]) is the energy flux going into the 
compartment i through the piston velocity fluctuations. This quantity has been evaluated 
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in 17| by means of a simplified kinetic theory in which both the velocity distribution of the 
piston and of the inelastic gases were approximated by Gaussians. Particularization of Eq. 
(66) in [lT I for the present case of elastic collisions between the particles and the piston gives 
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where 



mTp 
~MTi 



(8) 



(9) 



and Tp is the temperature parameter of the piston defined from the second moment , < 
V*{t) >, of its velocity distribution by M < V%{t) >= T P . Of course, no macroscopic or 
thermodynamic interpretation is assigned to this quantity, although it is expected to be 
a decreasing function of time. More precisely, the energy balance equation for the piston 
yields 

-2S(Q 1 + Q 2 ). (10) 



dt 

Equations ([TJ) , (J3J) , and f flUj) form a closed set of four equations for the unknown xp, Ti(t), 
T 2 (t), and T P (t). 

To simplify the calculations and allow for an analytical solution, attention will be re- 
stricted in the following to those cases where m/M is small, so that, assuming that Tp/Ti is 
of the order of unity (something to be checked a posteriori), Eq. (jSJ) can be approximated 

*- 2 (v) 1/2 T^ 1/2 <»> 

The above expression is consistent with the intuitively expected behavior. Energy flux goes 
from the system with the larger temperature parameter to the system with the lower one. 

III. THE SYMMETRIC STATE 



It is convenient to use a dimensionless time scale s defined by 

\ P (ty 1/2 



ds = (2n) 1/ V 



1/2 d-1 



m 



dt, 



(12) 



where n = 2N/ SL X is the average number density of the whole system. Then, Eqs. ([3]) and 
(05} read 



dJ\ 
ds 



-CLiTi 



2Nd 



{Ti — T P ), 



(13) 



i = 1, 2, and 

^ = |pi-IW + |(T,-T,). (14) 

Here, 

ai =(-) a, bi =^) b, (15) 
with a and 6 being constant parameters given by 

a = (2 + d)2 1 / 2 T(d/2) ' (16) 
_ SS'm 

= ,i/2 ff i-i M • ^7) 

The temperatures 7\ and T 2 are not independent since they are related by Eq. ([1]). Therefore, 
the three equations (fT3 l) -( fl5l) can be replaced by the set of linear equations 



is V ai + ^J T1 + ^ TP ' <18) 
2 dTi ( b\ \ m b\ 
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ds 2 y r ' 2p 
where the relative density of the two compartments 



P-^i, (21) 
n 2 

has been introduced. To search for solutions of the above equations write 

H (s) = ie Xs , T P {s) = 6 P e Xs . (22) 

Then it follows that for the existence of solutions different from the trivial one 9\ = Op = 
it must be 

a ' + 2^ + A =( a ' + 2^ + A ")" 2 < 23 > 

and 

2 ("i + 2^ + A ) _ bx (24) 



& i(!+P) 2Nd 



bi(l+P) 



+ A 



2p 

Consider first the special case p = 1, i.e. the symmetric state with n\ = n 2 = n and 
xp = L x /2. Because of Eq. ([T]) it is also 7i(£) = T 2 (t) = T(t) and, because of Eq. (jTTI) . 
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Qi — Q2 = Q- For this value of p, Eq. (1231) becomes an identity, and only the condition 
given by Eq. (I24jl remains. The latter takes the form 

X2+ { a + b+ 2Wd) X + ab = °' (25) 

and it has the solutions 

A± _c ± ( £i ^ 1 V !i (26) 

^° + 6 ( 1 + ^)' <27> 
It is easy to verify that A_ < A+ < and, therefore, for large enough times the temperatures 

of the granular gases and the piston in the symmetric state have the form 

r«( s ) = 9^e x(a)s , T p s \s) = ^ s) e A(s)s , (28) 
respectively, with A^ = A + . The coefficients 6^ and 9^p are related through 

that follows from Eq. (120]) after particularizing it for the symmetric state. Using Eq. (|26|) it 
can be checked that |A^| < b, as required by consistency, since otherwise the theory would 
predict negative values for one of the temperature parameters. In addition, Eqs. ( 128]) and 
(125)1 imply that 

T p s \t)>T {s \t), (30) 

i.e. the temperature of the granular gases is smaller than the temperature parameter of the 
piston at the same time. This reflects that the physical mechanisms for which the piston 
cools is the energy dissipation in the gas collisions, so that Tp is driven by T^ S K 

Equations (128]) indicate that the temperature parameters of both the piston and the 
granular gases cool at the same rate, 

^0) T / \ oip_ _ _As) T ( s ) 



dT^ s dTi s) 



V = -w T(s) ' ^ = <f T > (31) 

where 

< S ) -- A( "0 V>nff " (32) 

is the effective cooling rate of the gas taking into account the energy interchange with the 
piston. In the limit N — > 00 with all the other parameter remaining finite, it is A^ ~ —a 
and QJ ~ (, as expected. In the same limit, Eq. ( 1291) leads to 

(S) « (l " f ) 0p> (33) 
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showing that the symmetric state does not exist if a > b in this limit. But for M/m 3> 1, 
it is b <C 1. This means that for iV — > oo, M/m — > oo, the state with xp = L x /2 does not 
exist. The physical reason is that when the piston has a very large mass and the number of 
particles is very large, the piston is not able to cool down with the same rate as the granular 
gas. 



IV. THE ASYMMETRIC STATE 



To investigate the existence of states with a stationary position of the piston different 
from the middle of the system, Eqs. (123]) and (|24|) have to be solved. For p =fi 1, eq. ( 123]) 
yields 

In the above expression p has been eliminated in favor of the number densities to emphasize 
the invariance with respect to the interchange n\ -h- ri2, as required by symmetry. Therefore, 
it is 

if = e[ s) e x(a)s , T 2 (a) (s) = 0£V (a \ (35) 
T p s) = ^ } e A(<l)s , (36) 



with 
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(a) ^ ^ 

^ = 1 1 + ^J7T7TT- (38) 

The last relation follows directly from Eq. ([18]) . Still, it must be required that Eq. (124]) 
be fulfilled but, before doing that, it is important to realize that Eqs. (135]) and (136]) lead 
to the result that the cooling rate of both granular gases and also the piston are again the 
same in the asymmetric state(s), namely 

^_ = _ c W T W (t)j (39) 



dT[ a) 



= _ c W T W (t)> (40) 

cg^-A^f^y 72 ^- 1 . (4i) 
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Substitution of Eq. (1341) into Eq. (|24|) . after some algebra leads to 



z 1 - Ndz + Nd = 0, 



with 



(IN ad 



P 



p 2 + p + 1 



For large values of N, the two solutions of Eq. (1421) read 



Z\ « 1, z 2 ~ Nd. 



(42) 



(43) 



(44) 



The next task is to invert Eq. (14*3]) to get the values of the density ratio p. Consider first 
the root z\ — \. Because of Eq. f|38|) it is 

1/2 



or, in terms of the actual temperature parameters, 



Tf(t) = T[ a ' '(t)T 2 w \t) 



(o)/.^(a). 



1 1/2 



For z — 1, Eq. ff43|) becomes 



j(a)2 (a) + 1 = 



whose solution is 
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- 1 



(45) 



(46) 



(47) 



(48) 



(The other solution of the second degree equation is the inverse of the above and corresponds 
to the interchange of the compartments 1 and 2). Existence of a physical solution requires 
that 



Nad 



> 1. 



Using Eqs. ( IT6|) and ( TT7|) . this condition can be expressed as 



M > M r 



where 



M r EE 



(ci + 2)2 1 / 2 r (d/2)Sm 



(49) 



(50) 



(51) 



d7r d/2 a d-l N £* 

can be understood as a critical value of the mass of the piston below which the asymmetric 
state with z — 1 is not possible. Moreover, the larger M/M c the stronger the asymmetry of 
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the state as measured by the position of the piston xp. In the limit M/M c — >■ oo , either 
xp — > ot xp — > L x . This is just the limiting situation discussed in ref. p. 

It is relevant to see what happens at the bifurcation point, i.e. for M — > M c or, more 
precisely, for M - M c -»■ + . In this limit p -»• 1, and Eqs. (j3"5) - ([3gj) yield 

T[ a \t) =T^ a \t) =T p a \t). (52) 

Moreover, Eq. ( 134|) gives — >■ —a. This value agrees with the relaxation rate A'- 5 ) of the 
symmetric state, given in Eq. ( 12 6 p in the limit Nd 3> 1, indicating the continuity of the 
cooling rate through the transition. 

Still remains to be studied the solution Z2 ~ Nd of Eq. fj42l . It is easily seen that this 
value leads to the relation 

T P (t)^' = Nd [T/ a)/ (t)T 2 (a)/ (t)l V2 , (53) 

following from Eq. ( )38|) . It follows that this solution has to be discarded from our analysis 
since for large N the temperature parameter of the piston is much larger that the temperature 
of the granular gases to both sides of it, something that contradicts the assumption made to 
derive the approximate expression for Qi given in Eq. (fTTj) and used throughout the paper. 

V. MOLECULAR DYNAMICS SIMULATION RESULTS 

To check the above theoretical predictions, molecular dynamics (MD) simulations with an 



event driven algorithm [18l . Il9| have been performed. The system considered is an ensemble 
of 2N inelastic disks enclosed in a two-dimensional rectangular box of sides L x and L y . 
The latter corresponds to the general transversal section S used in the previous sections. 
Collisions of the particles with the walls of the container and also with the piston are elastic. 
In all the simulations, the piston was initially placed in the middle of the container and 
the particles were evenly and uniformly distributed between the two compartments. The 
initial number density was always quite low. Moreover, the initial velocity distributions of 
the disks were Gaussian and with the same temperature at both sides of the piston. On the 
other hand, the latter was placed with a vanishing initial velocity. 

There is a limitation in the systems that can be simulated as a consequence of the 
instability of the HCS. As the value of the coefficient of normal restitution a decreases, the 
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20| . that means that 



critical size for which the HCS becomes unstable also decreases [13 
for a given density the number of particles becomes smaller. In the low density limit, this 
leads very soon to a very reduced number of particles, so that boundary effects becomes 
very relevant, specially when dealing with a system having walls (as opposite to periodic 
boundary conditions), as it happens in the present case. As a consequence, it turns out that 
the particular states addressed here are not reached for values of a roughly smaller than 
0.85. 

In some cases, it was observed that the piston stays oscillating around the initial position, 
xp = L x /2. These oscillations are easily understood as being induced by the pressure 
fluctuations of the granular gas. This behavior persists in time, so that the initial symmetry 
of the system is conserved. It is clear that these situations correspond to the symmetric 
state discussed in Sec. IHIl On the other hand, in other cases the initial oscillations increase 
in time and after a while they combine with a net motion of the piston moving away from 
the middle of the system. Eventually, the net motion ceases and a steady average position 
is observed. The system reaches one of the asymmetric states analyzed in Sec. IIVI Figure [2] 
shows a typical example of the observed behavior when the initial symmetry is broken. The 
parameter values in this case are a = 0.98, 2N = 200, L x = 2L y = 100a, and M = 25m. 
Then, the initial number density is n\ = ri2 = 0.04cr -2 . Note that the time scale r used in 
the figure is the accumulated number of collisions per particle. In the particular simulation 
reported, the piston moves towards the right compartment, reaching a stationary average 
position at xp ~ 0.66L X .. In different simulation realizations with the same values of the 
parameters, the final position of the piston occurs equally often to each side of the middle, 
but always at the same average distance of it. Moreover, it is worth to mention that in some 
realizations spontaneous transitions between the two average positions at both sides of the 
middle were observed. 

When the average position of the piston reaches a steady value, also the average density 
profile becomes time independent. It is given in Fig. [3] for the same values of the parameter 
of the system as in Fig. [2j The plotted curve is an average on time once in the steady con- 
figuration and also over 350 simulation trajectories. It is observed that outside a boundary 
layer at both sides of the piston, the density can be considered as uniform and different in 
each of the two compartments, in agreement with the assumption made in the theoretical 
model developed in Sec. [Til The density boundary layer is mainly due to the oscillations of 
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FIG. 2: Time evolution of the position of the piston for a particular simulation trajectory for a 
system with 2N = 200, a = 0.98, L x = 2L y = IOOct, and M = 25m. Time r is measured in 
accumulated number of collisions per particle. The dashed line is a guide for the eye. 
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FIG. 3: Steady density profile exhibited by the the system with the same parameters as in Fig. (2) 
The solid vertical line indicates the steady average position of the piston. 

the piston. 

Another strong prediction of the theory is Eq. ( I46I) . relating the steady temperature of 
the piston and those of the inelastic gases in both compartments. The MD simulations show 
that the relation is verified in the range of parameters in which the asymmetric state has 
been observed. In Fig. HJ this is illustrated for the same values of the parameters as in Figs. 
[2] and [3j To put the comparison in a proper context, it must be taken into account that 
along every trajectory the piston jumps several times from the average position at one side 
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FIG. 4: Temperature of the piston T p divided by the square root of the product of the temperatures 
of the granular gases in each of the compartments, T\ and T2, respectively, for the same system 
as in Figs. [2] and [3l The data have been averaged over 2000 trajectories. Time is measured in 
accumulated number of collisions per particles. 

of the system to the symmetric one, through states which are not described by the theory 
developed here. These configurations have not be removed when computing the results 
reported in the figure. It is worth to notice that Eq. (l4"6"j) is consistent with the result that 
the cooling rates of the gases and the piston are the same. 

Therefore, the simulations clearly show the existence of both the symmetric and the 
asymmetric states. In the latter, the steady average position of the piston can take any 
value inside the system, depending on the values of the parameters. In order to carry out 
a more quantitative check of the theory, the average position of the piston xp has been 
measured as a function of the ratio M/M c , where M c is the parameter denned in Eq. (I5ip . 
The results are plotted in Figs. [5] and |6] in terms of the asymmetry parameter 



\2x t 



Il-P 2 



(54) 



L x 1 + p 2 

defined in the interval < e < 1. The prediction for this quantity obtained here is 

^ = 1 and e = for M < M c (55) 



P = P 



P = P 

(a) 



M 



M 
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1/2 



for M > M r . 



(56) 



The observed agreement between theory and simulations can be considered satisfactory. 
In particular, the bifurcation point seems to be quite accurately predicted by the theory. 
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FIG. 5: Bifurcation diagram showing the asymmetry of the number density as a function of the 
dimensionless control parameter M/M c , where M is the mass of the piston and M c is given by Eq. 
(|5ip . The symbols are from molecular dynamics simulations and the solid line is the theoretical 
prediction derived in the paper, Eqs. ([55]) and {(56]) . All the simulation results reported in this 
figure have been obtained with L x = 2L y = 100<7 and 2N = 200, so that na 2 = 0.04. Different 
values of a have been used as indicated in the insert. For each value, the mass ratio M/m has 
been varied in order to change the value of M/M c . 

Upon evaluating the results, it must be taken into account that when the asymmetry of the 
state increases the density in one of the compartments also increases, and that the theory 
developed here is based on the assumption that the granular gas in both compartments can 
be treated in the very dilute limit approximation. 



VI. CONCLUDING REMARKS 



The spontaneous symmetry breaking investigated here has some similarities with phe- 
nomena of phase separation occurring in vibrofluidized granular materials, both in presence 
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22] and absence 23] of gravity. There are also significant differences. A relevant one 



is that the system studied in this paper is isolated and has no stationary states, as a con- 
sequence of the energy dissipation in collisions. The states addressed are characterized by 
a stationary position of the movable piston dividing the system into two parts. A main 
result obtained is that the conditions determining these states do not require the (granular) 
temperature to be the same at both sides of the piston, as it is the case in normal, elastic 
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FIG. 6: The same as in Fig. \5\ with the only difference that in the simulations reported here it is 
L x = ALy = 200cr and, consequently, no 2 = 0.02. 

fluids. Instead, the characteristic relaxation times of the energy must be the same in both 
compartments, and agree with the one of the piston. These times can be interpreted as the 
effective cooling rates of the granular gases. Moreover, the gases tend to be homogenous at 
both sides of the piston, although presenting boundary layers next to it. 

There is an interesting physical picture emerging from the above results. It is well known 
that when two molecular systems are put into contact through a movable thermal wall, 
keeping the whole ensemble isolated, the final state is of equilibrium, implying that each of 
the two macroscopic subsystems is also at equilibrium, having both the same pressure and 
temperature. If the same experiment is carried out with granular gases, the results being 
reported indicate that the whole system tends to a cooling state characterized by a unique 
time-dependent temperature parameter T(t) and being homogeneous inside each compart- 
ments. To be more precise, what is meant is that at a microscopic level the probability 
distribution function of the system has the form 

/<™ = m«)]-"™' -Mi) }•!))• (57) 

where T denotes a point of the phase space of the system, T{ and (i = 1, . . . ,2N) are 
the position and velocity of particle i, X and V x are the position and velocity of the piston, 

1 /2 

and v(t) = (2T(£)/m) 1 . The peculiarity of this state is that all its time dependence 
occurs through the parameter T(t). The function /* includes the constraint that some given 
particles are located at one side of the piston and the remaining particles at the other one. 
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Temperature parameters Tj(t), i = 1,2, for each compartment can be defined in the usual 
way from the second velocity moment of its particles. Also a temperature parameter Tp(t) 
for the piston can be defined in a similar way. The scaling form of Eq. (1571) implies that 

dlnTi(t) d\nT 2 (t) d\nT P (t) 



dt dt dt 



(58) 



i.e., the cooling rates of both compartments and also of the piston are the same. This 
is the result found in this paper. Trivially, Eq. ( |58l) is verified if all the temperatures 
are the same, i.e. T\(t) = T 2 {t) = Tp(t). This is what happens in equilibrium molecular 
systems, where velocity correlations between the three macroscopic systems are negligible 
and the equilibrium distribution for the velocity factorizes. On the other hand, when dealing 
with granular gases, velocity correlations between the two compartments and the piston are 
relevant, the velocity distribution of the whole system does not factorize and the three 
temperature parameters differ. 

The above feature is related with the non-equipartition of kinetic energy in mixtures of 
granular gases. The granular temperature of the components of a mixture defined from the 



average kinetic energy of each species is not the same [2J, |25|. The homogeneous cooling 
state of the mixture has the property that the cooling rates for all the partial temperatures 
are the same. This condition determines all the temperatures of the components in terms 
of a unique temperature parameter. Actually, a similar result holds in vibrated granular 
systems, in the sense that only a temperature is needed for a macroscopic description of 
mixtures of granular gases [26j. Here, the result is extended to spatially separated granular 
gases interacting through a movable piston. 

A consequence of the cooling rate being the same as compared with the temperatures 
being the same, is the existence of the spontaneous symmetry breaking discussed in this 
paper. The phenomenon can be easily understood in terms of rather simple kinetic theory 
arguments, whose theoretical predictions are in good agreement with molecular dynamics 
simulation results. An interesting issue to be addressed in the near future is why the asym- 
metric state is more stable than the symmetric one when both states exist. Perhaps this is 
related with some extremal principle that shed light on the macroscopic physics of granular 
fluids. 
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